Phase of the Fermion Determinant at Nonzero Chemical Potential 
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We show that in the microscopic domain of QCD (also known as the e-domain) at nonzero 
chemical potential the average phase factor of the fermion determinant is nonzero for n < m7r/2 
and is exponentially suppressed for larger values of the chemical potential. This follows from the 
chiral Lagrangian that describes the low-energy limit of the expectation value of the phase factor. 
Explicit expressions for the average phase factor are derived using a random matrix formulation of 
the zero momentum limit of this chiral Lagrangian. 
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Introduction. During the past decade, a great deal 
of progress has been made in understanding the phase 
diagram of the QCD partition function in the chemi- 
cal potential - temperature plane. Although early an- 
alytical arguments [jj] clarified the nature of the chiral 
phase transition along the temperature axis, a detailed 
quantitative understanding could only be achieved by 
means of lattice QCD simulations (see j^j for a review). 
The situation at nonzero chemical potential is much less 
clear. Although perturbative arguments, model calcu- 
lations and phenomenological arguments seem to give 
a consistent picture Q, first principle quantitative in- 
formation is lacking. The main reason is that QCD at 
nonzero chemical potential cannot be simulated reliably 
in much of the chemical potential - temperature plane be- 
cause the fermion determinant is complex (the sign prob- 
lem). Progress has been made around the critical tem- 
perature and small chemical potentials, where different 
lattice QCD approaches seem to converge 
Throughout this letter, chemical potential is short for 
quark chemical potential and will be denoted by /x. 

The question we wish to address in this letter is 
whether there exists a parameter domain for which the 
phase of the fermion determinant is manageable. Since 
there is no sign problem for /i = 0, we expect that 
for sufficiently small nonzero chemical potential lattice 
QCD simulations are possible. The standard argument is 
that the number of gauge field configurations required for 
a converged calculation diverges exponentially with the 
volume as exp(FAF), where VAF is the difference of 
the free energy of typical gauge field configurations gen- 
erated by the Monte Carlo algorithm and the converged 
free energy. One limit in which the sign problem remains 
manageable is when VAF remains finite in the thermo- 
dynamic limit. We expect that for sufficiently small fi 
the free energy difference behaves as fi^V. For the sign 
problem to be manageable the chemical potential then 
has to scale as 1/VV in the thermodynamic limit. Such 
type of limit is well known - it is the microscopic limit or 
the e-limit of the QCD partition function. In this limit 
the mass and chemical potential dependence of the QCD 
partition function is determined by the zero momentum 



modes of the Goldstones associated with the spontaneous 
breaking of chiral symmetry. 

One observable that directly tests the severity of the 
sign problem is the expectation value of the phase factor 
of the fermion determinant. This phase factor can be 
expressed as the ratio of the fermion determinant and its 
complex conjugate 
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Its expectation value is a QCD-like partition function 
with a low energy limit that is determined along well- 
established rules by chiral symmetry and gauge invari- 
ancc 9, 10]. In this letter we analyze the average phase 
factor in the microscopic domain of QCD JJ^ 
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In this domain the Compton wave length of the Gold- 
stone modes is much larger than the linear size of the box, 
so that the chiral Lagrangian can be truncated to its zero 
momentum sector. In the thermodynamic limit, simple 
expressions can then be obtained using mean field argu- 
ments. At finite volume, the calculations are much more 
complicated, but we can exploit the equivalence with 
random matrix theory [l^ [l^ . where recent progress 
El 111 in III III Eg makes it possible to derive ex- 
act results in the microscopic domain. Several cases will 
be discussed: the quenched limit, the phase quenched 
limit, and QCD with dynamical flavors. In all cases we 
will flnd that the average phase factor is nonzero even 
for large fi'^V provided that 2fi < TOtt. For 2^ > rriT^ the 
average phase is exponentially suppressed with ^^V. 

Although QCD at nonzero baryon chemical potential 
has a sign problem, this is not the case for QCD with 
two colors, QCD with gauge fields in the adjoint rep- 
resentation and the phase quenched partition function. 
The chemical potential and mass dependence of these 
partition functions has been analyzed in great detail by 
means of chiral lagrangians or random matrix theory 
llElM m 0, Q 123 as well as on the lattice 
123, 12911 . The success of these calculations suggests 
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that equally impressive lattice QCD results can be ob- 
tained for the average phase factor. We hope that the 
results presented in this letter will encourage such calcu- 
lations. 

The approach introduced in this letter is directly appli- 
cable to QCD at nonzero 0-angle. Fermion sign problems 
also appear in other interesting physical systems • It 
would be worthwhile to analyze them along the lines pro- 
posed in this letter. 

General arguments. In this section we will evaluate the 
/i-dependence of the average phase factor in the mean 
field limit. Below we will confirm these results from the 
asymptotic limit of the exact expressions for the micro- 
scopic domain. In this domain, it is natural to work at 
fixed topology instead of fixed 0-angle. The results pre- 
sented in this section are for the thermodynamic limit 
and do not depend on the topological charge. 

The vacuum energy density does not depend on the 
chemical potential in a phase that is not sensitive to the 
boundaries. This is the case in the normal phase where 
the chemical potential is below the mass of the light- 
est particle with the corresponding charge. For larger fi 
there is a net particle flux in the time direction of the Eu- 
clidean torus and the free energy depends on the chemical 
potential ;31]. Although in the normal phase the free en- 
ergy is /X- independent, the excitations of the vacuum are 
not. For a chiral Lagrangian the masses of the Goldstone 
modes for fi < are given by ^3 



M{p) — m^r — b^, 
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where b is the charge of the particles corresponding to fi. 
In the zero momentum sector, the thermodynamic limit 
of the partition function is therefore given by 
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where the Jacobian, J, is from the measure of the Gold- 
stone manifold and F is the vacuum energy density, both 
evaluated at the saddle point. The prefactor gives a 1/V 
correction to the free energy density. The prefactor is im- 
portant if we consider the expectation value of the phase 
factor of the quark determinant which is given by the 
ratio of two partition functions 



(e iNf = 



ZNf 



(5) 



They are defined by ((••■) refers to quenched averaging) 



/ det(Z) + -I- m) jy 
' \ det(L'T -I- ^-)/Q -I- m) 



IJ-lo + m)) ,(6) 



The partition function ((BJ contains |33 Nf + 1 fermionic 
quarks and one conjugate bosonic quark. Assuming max- 
imum spontaneous breaking of the axial flavor symme- 
try, this results in Nf -\- 1 charged fermionic Goldstones 
composed of a fermionic quark and a conjugate bosonic 
anti-quark as well as an equal number of anti-particles 
with the opposite charge. In addition, for topological 
charge zero, we have the usual {Nf + Vf neutral bosonic 
Goldstones and one neutral Goldstone made out of two 
bosonic quarks. The partition function in the denom- 
inator contains N'j neutral Goldstones. In the normal 
phase, the saddle point of the static part of the effective 
Lagrangian is /^-independent and neither the Jacobian 
nor the free energy do depend on /x. Using |0J, the aver- 
age phase factor for ^ < TOtt/Z is given by 
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We emphasize that the free energies of Z7Vj>+i|i» and Zffj 
cancel. Hence for /x < TO,r/2 the sign problem is not 
exponentially hard in the microscopic domain 1^)1. 

This is not the case for /i > m7r/2, where the free 
energy of Z^^^i^i. is /x-dependent exactly as in other 
theories with charged Goldstone particles [lOj . This leads 
to an exponential suppression of the average phase factor 
{Ft: is the pion decay constant) 



(9) 



for /i > m-^/2. In addition, the Goldstinos with mass 
TOtt — 2/i become exactly massless for /x > 771^/2 so 
that the leading contributions to the prefactor cancel. 
We conclude that the sign problem is not tractable for 
large ^J.'^F^V and fi > 

Microscopic result. As we have seen in Eq. (jSl, the 
expectation value of the phase factor is given by the ratio 
of two partition functions. We now calculate them in the 
microscopic limit where the scaling variables 



and fi = ^iF-^Vv 



(10) 



are kept fixed for V oo. In this limit, the QCD parti- 
tion function is equivalent to the large N limit of a ran- 
dom matrix theory of 2N x 2N matrices with the same 
global symmetries and transformation properties 
This allows us to perform the calculations using recent 
developments in the method of orthogonal polynomials 
llTj J9(,,22j|. Starting from a general expression in [lH, 
it can be shown that the microscopic limit of the parti- 
tion functions in Q can be expressed in terms of modified 
Bessel functions and their Cauchy transforms. For zero 
topological charge we obtain (with 6,% = md/dm) 



and 



ZMf = {dei^f{D + ^iJo + m)). 



(7) 
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where llll 



. ^-Nf(Nf-i) det[4+'/o(m)]fc,,^o...iv,-i. (12) 



The Cauchy transforms X^'''>{m; jl) are defined by 



d^z 



^(z, z*, /o(z* 
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where i(j(z,z*;/i) is the weight function of the random 
matrix model in [l^ . The expressions for the Cauchy 
transform can be rewritten as a one-dimensional in- 
tegral following the approach of Next we give ex- 
plicit results for the thermodynamic limit of Hll|l which is 
obtained from the saddle point approximation for jl —t oo 
and m — > c» at fixed fj? /m. 

In the quenched case {Nf = 0) a saddle-point approx- 
imation of Hll|l gives 

(e2'«)Ar,^o = (l-^)e", 2^1 <m^. (14) 



This result agrees with the mean field arguments given 
above. For /i > m7r/2 the result is exponentially sup- 
pressed exactly as in (O and with a prefactor that can- 
cels to leading order in \/V . Notice that the phase factor 
is only exponentially suppressed for /i > 

For Nf = 1 the thermodynamic limit of the exact mi- 
croscopic result (|llf) is given by 



(e )Nf=i 



(1 - Kf 

mi 



2/i < m^, 



(15) 



in agreement with the mean field arguments given above. 
For jl > at finite volume the result is given by © 

with a prefactor that cancels to leading order in 1/V. 

By now it should be clear that the thermodynamic 
limit of the microscopic result 1)11(1 reproduces the general 
formula ||SJ| for all values of Nf. As further illustration 
we plot in Fig. 1 the average phase factor for Nf = 2. 
The dashed curve represents the result of Eq. ifTT)) for 
mVT, = 4 and the full curve is its limit for mVY, oo 
at fixed ///mTT- We observe a rapid convergence to the 
thermodynamic limit especially at small 

Finally, we calculate the average phase factor for the 
phase quenched theory where the phase factor of the dy- 
namical fermions is ignored. For two fiavors it can be 
expressed as 



■lie 



(|det(i? + /i7o-Hm)|2) 



(16) 
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FIG. 1: The average phase factor for two flavor QCD in the 
e-regime as a function of the chemical potential for m'EV — 4 
(dashed curve) and m'EV oo (full curve). 



The microscopic limit of both partition functions is well- 
known [llj, |23| , resulting in the average phase factor 
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In the thermodynamic limit obtained by making a saddle 
point approximation of the integrals in H17|l we find the 
same result as in the quenched case (|14f) . For 2/z > 
the result is once more given by (^J but with a different 
prefactor than in the previous two cases. 

Lattice Simulations. Several lattice simulations have 
studied the fluctuations of the phase of the fermion de- 
terminant. Both {cos 9) ,.32j and {6'^) 'ss'l were consid- 
ered. In the variance of the phase of the fermion 
determinant was calculated using the Taylor expansion 
technique. To lowest order in an expansion in /i the av- 
erage squared phase is given by 



M^t^Mi^M2 log ^1+1* 



Imi=M2=0 ' 



(18) 



where 



n+i' 



(det(D -I- ^i7o -I- m) det(D'l' + /i27o + m)).(19) 



In the microscopic limit, this partition function is given 
by the the denominator of Eq. 1(17(1 with u — > (/ii — /i2)/2. 
In the thermodynamic limit we obtain {/j, < m,r/2) 



(20) 



The work of AUton et al. indeed suggests that 

(6*2) ~ f^'^/ml- but the prefactor appears to be several 
times larger than given by 1(20(1 . There can be several 
reason for this discrepancy. First, the calculations of All- 
ton et al. 1^ were performed close to the critical tem- 
perature whereas our results are for zero temperature. 
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In particular, susceptibilities are expected to be sensitive 
to the temperature. Second, our results have been de- 
rived for the e-domain of QCD whereas the pion mass 
in 1^ does not satisfy the condition Third, there 

could be significant ultra-violet contributions to the av- 
erage squaredphase. Although, it can be shown along 
the lines of [26j that in dimensional regularization the 
^—dependent terms do not introduce additional ultra- 
violet divergences, for a lattice regularization this is only 
the case after the necessary subtractions have been made. 
Ultra-violet contributions to the average phase factor are 
expected to behave as exp(— Fu>^A^) with w the width 
of the strip of eigenvalues and A an ultra-violet cut-off. 
Since w ^ fj? (see 0|) ultra-violet contributions are sup- 
pressed in the microscopic limit. 

Conclusions. We have shown that for sufficiently small 
fj, the expectation value of the phase factor of the quark 
determinant can be obtained from chiral perturbation 
theory. Explicit expressions have been obtained in the 
microscopic domain where /i ~ 1/VV and m ~ ^/V. Our 
results show a phase transition of (exp(2i0)) at fi — 
with (exp(2i6')) = beyond this point. In particular, 
this implies that there is no serious sign problem for /x ^ 
and a physical quark mass. However, the sign 
problem is severe in the chiral limit for any nonzero /i 
where it is essential for the discontinuity of the chiral 
condensate [s^ . 

Our results have been derived for zero temperature. 
From the temperature dependence of the grand poten- 
tial we expect that phase fluctuations initially increase 
with temperature. A deeper understanding of the sign 
problem could be obtained by extending the current lat- 
tice simulations to lower temperatures and quark masses. 
We believe that the confirmation of the analytical results 
for the phase fluctuations predicted in this letter will be 
an important step forward toward a first principles under- 
standing of the QCD phase diagram at nonzero chemical 
potential. 
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